Cluster Spin Glass Distribution Functions in La 2 x Sr x Cu04 
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Signatures of the cluster spin glass have been found in a va- 
riety of experiments, with an effective onset temperature T on 
that is frequency dependent. We reanalyze the experimental 
results and find that they are characterized by a distribution 
of activation energies, with a nonzero glass transition tem- 
perature T g (x) < T on . While the distribution of activation 
energies is the same, the distribution of weights depends on 
the process. Remarkably, the weights are essentially doping 
independent. 



I. INTRODUCTION 

Hole doping in the cuprates is now widely believed to 
lead to strongly inhomogeneous phases, with the holes 
being ejected from the antiferromagnetic domains onto 
antiphase boundaries. Ordering of these charged stripe 
walls leads to a form of domain phase, which is typically 
dynamic. In the very low doping limit, the inhomogeneity 
is manifested in the form of transverse spin freezing in the 
antiferromagnetic phasc,El and, when the Nael tempera- 
ture Tn — * 0, a^a 'cluster spin glass' (CSG).El Wakimoto, 
et al. (WUEY)El have correlated a number of experimen- 
tal observations of this CSG in I^-ajSr^CuO^ (LSCO), 
and raise the issue of whether there is a real finite temper- 
ature glass transition, or just a continuous relaxational 
slowing down. On the one hand, all of the experiments 
find the transition at a different temperature, depend- 
ing on the time scale of the measurement; on the other 
hand, the susceptibility shows striking scaling behavior, 
both with temperature and applied field. 

Here, we analyze these data in more detail, combining 
them with additional measurements at intermediate fre- 
quencies. In all measurements, we find that the CSG 
freezing turns on at an onset temperature T on which 
scales with the logarithm of the measurement frequency, 
consistent with a relaxational form, 
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with Vogel-FulcherD relaxation 

r = r e E ^ T - T ^ 
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(here T on corresponds to the temperature below which 
t -1 becomes slower than the measuring angular fre- 
quency lj). In order to explain the observed scaling be- 
havior, we find that the CSG must have a broad distri- 
bution of activation energies, with a well defined charac- 
teristic energy scale E m , and that different experiments 
probe different aspects of this distribution. 



II. T G AND ACTIVATION ENERGY 
DISTRIBUTION FUNCTIONS 

The glass transition temperature T g is typically found 
by looking for deviations from activated behavior: cur- 
vature in a plot of ln(oj) vs. 1/T (here, u> is a char- 
acteristic frequency of a given experiment and T is the 
temperature). The presence of a broad distribution of 
activation energies makes it difficult to extract a unique 
value of T g . For the CSG, there is an additional com- 
plicating factor: the characteristic frequencies are often 
not well defined. p. For instance, in the neutron diffrac- 
tion experiments^ the effective frequency is taken as the 
energy resolution 0.25meV: the diffraction peak is elas- 
tic, or static, on that energy scale. On the other hand, 
the magnetic susceptibility was measured by a 'static^ 
SQUID technique since a typical scan takes about 10 s,cl 
we assume an effective frequency of ~ 1/20 Hz. Clearly, 
improved determination of these frequencies will lead to 
better estimates for T g - particularly at the lowest fre- 
quencies. Nevertheless, the variation with frequency is 
so striking that the qualitative features should be un- 
changed by these refinements. 

The simplest distribution comes from measuring the 
magnetic susceptibility \oiT). WUEY found that the 
in-plane xo could be fit to a simple Curie law with small 
Curie constant for temperatures between 10 K and 70 K, 
while at lower temperatures the Curie constant became 
temperature dependent as clusters began to freeze out. 
To analyze these results, we assume that the susceptibil- 
ity has a relaxational component 



Xo (T) = xc(T)(l -g)= dEw susc (E) 
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1 + lu 2 t(E) 2 ' 
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with xc the Curie susceptibility and w SU sc (E) the suscep- 
tibility distribution function. This distribution function 
has a simple interpretation in the standard model of 
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The above analysis can be extended to other exper- 
imental measurements. There are two changes. First, 
the susceptibility is proportional to the product of a 
coupling and a distribution function, xc(T) w susc (E); 
for other properties, this can be written Ai(T)wi(E), 
where Ai represents the intrinsic temperature depen- 
dence of the process. In these calculations, we assume 
that A ane i ast — 1/T (Ref. |g) while A NQR — const. Sec- 
ondly, since most of the other probes are sampling a dis- 
sipation, there is an extra factor of lot in the integrand, 
as in Eq.|l|. This factor means that the integrand at any 
T will be strongly peaked at E(T). 



FIG. 1. Distributions of excitation energies w S usc (solid 
line), w ane i a3t (dotted line), and w NQR (dashed line); 
dot-dashed line = Gaussian approximation to wnqr- Insert: 
resulting q(T), Eqj|. 

the susceptibility. It is assumed^ that there is a distri- 
bution of clusters, with their spins locked in an antifer- 
romagnetic pattern, so that each cluster with an even 
(odd) number of spins has net spin zero (1/2). The 
observed Curie law is attributed to the response of the 
odd-spin clusters, thereby naturally explaining the small 
Curie constant. Since the Curie law is obeyed between 
10-70 K, the clusters must have a well defined size dis- 
tribution, which begins to freeze out at lower temper- 
atures; w susc (E) is then the distribution of the num- 
ber of clusters with a given activation energy, indepen- 
dent of the cluster size. If this distribution is normalized 
(/ dEw susc (E) = 1), then 
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The strong temperature dependence of r, Eq.^J, means 
that Eq. |j can be simplified. For each temperature, 
there is a corresponding activation energy for which 
ujt{E{T)) = 1, namely E(T) = -Tln(u>T ). Then the 
effective Edwards- AndersonB order parameter q may be 
approximated 



(i / dEw susc (E). 

I E(T) 



(5) 



The observed^ scaling q ~ (T on — T) 13 with /3 ~ 1 then 
suggests that w susc {E) — const for E < E m , with E m — 
E{T on ). Thus the appearance of T on is interpreted as the 
presence of a sharp cutoff E m , 

T on = -E m /ln(wT ). (6) 

The (flat) distribution of w susc is shown in Figure [l]; the 
insert to the figure shows that the above analysis is well 
satisfied by direct numerical integration of Eq. ||. We 
note that Eq. |^ has been applieda for describing suscep- 
tibility data in the spin glass state (T < T g ); however, 
the above analysis depends only on the form of the dis- 
tribution, and holds even if T g = 0. 




FIG. 2. (a) Scaling of effective transition temperature (cor- 
responding to largest activation energy) with frequency. The 
symbols represent the measurement techniques: squares = 
neutron diffraction; circles = NQR; filled triangles = /iSR; 
open, inverted triangles = anelastic relaxation (1st mode - 
open triangle = 5th mode); and diamonds = susceptibility. 
The different sets connected by dashed lines refer to individ- 
ual doping densities, from left to right: x — 0.02, 0.022, 0.025, 
0.03, 0.04, 0.05, and 0.06. Solid lines = fits. Insert (b): corre- 
sponding T g (diamonds) and E m (x)/8 (circles), with dotted 
line = fit. T g is multiplied by 8 for comparison with E m . 

Figure |^a summarizes the doping dependence of T on 
vs. measurement frequency. The data are, from highest 
frequency: squares = neutron scattering determination 
of the onset of (diagonal) stripe orderQ (ui =j— 0.25 meV 
1 s~ 1 ); circles = NQR measurementstil (w/27r = 
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19 MHz); filled triangles = fiSaB (u = 10 6 s" 1 ); open 
triangles = anelastic relasationliil (uj/2ir ~ 2 kHz); and 
diamonds = susceptibilitya (u/2tt ~ 1/20 Hz). Solid lines 
are fits assuming a finite glass transition temperature T s ; 
Figure |^b shows the parameters T g (E m ) assumed in the 
fits. [In order to plot Fig. ||, it was sometimes necessary to 
use the smoothed curves of Ref. || rather than the actual 
data; this would only be a problem in the doping range 
x < 0.025, where some curves are extrapolated. Also, 
two fits are shown for the x = 0.06 data: the straight 
line corresponding to the larger E m and smaller T g (= 
0).] While there is considerable scatter in the data, it 
appears that To is essentially doping independent, while 
E m decreases with doping, Fig. ||b. For most dopings, 
the approximate relation E rn = &k^,T g is satisfied. The 
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value of t 1 = 5 x 10 13 s 1 = 33 nielLis close to the 
value 22 meV estimated by Julien et al£3 from the peak 
in the NQR Tf 1 . The dotted line in Fig. |b is a fit to the 
form E m = E /(x - x ), with E = 0.564 K, x = 0.011. 
The general trend E m ~ 1/x has been noted earlier .til 

We have included in Fig. || some anelastic and NQR 
data for a sample with x ~ 0.02. Uncertainty in the pre- 
cise doping makes it unclear whether the actual doping 
is x ~ 0.017, with spin freezing transition Tf ~ 13.9 K or 
x > 0.02, with T on ~ 10 K. We find that the assumption 
of an effective density x e ff = 0.021 yields remarkable 
agreement for both T on {u), but also for the weight, dis- 
tributions for both anelastic relaxation and NQRE3, as 
illustrated in-Figs. [| and || below. Indeed, results by 
Curro, et al£3 show that La NMR weight distributions 
are very similar in the spin freezing and CSG regimes. 
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FIG. 3. Elastic energy loss coefficient, Q -1 (T), for a x — 
0.03 sample, by exciting the lowest three (odd) flexural modes: 
squares = 1st (ui/2ti = 1.7 kHz), circles = 3rd (9.3 kHz), 
triangles = 5th (23 kHz). Solid lines = fit to the first mode, 
with distribution w ane i a3 t shown in Fig. 1, scaled by ui for 
the higher modes. 

Figure || illustrates the fit of the elastic energy loss co- 
efficient Q^ 1 for the three lowest odd modes of an LSCO 
sample at x = 0.03, while the corresponding distribution 
function is shown in Figure [j]. Approximately the same 
distribution works for the x ~ 0.02 and the x = 0.06 
samples (Fig. |^). Figure || shows the corresponding fits 
for the NQR relaxation rate; the distribution function 
(shown scaled in Fig. [l]) was fit for the x ~ 0.02 sam- 
ple, then scaled for the others. Note that for the anelas- 
tic response the fit is worse for the higher harmonics, 
and for x ~ 0.02 anelastic data there is an additional, 
lower-tcmpcpature process with a very different doping 
dependence.tH Also there appears to be a weak extra fea- 
ture in the NQR (~ 0.02) data near 20 K. It should be 
cautioned that the inversion of the data to obtain the 
distributions, Fig. [|, is not unique: the lorentzians are 
sufficiently broad that several distributions could lead to 



the same final spectra. The distributions chosen in Fig. [T] 
were chosen to be relatively smooth. [Note that the low- 
energy cutoffs are exaggerated, since data at sufficiently 
low temperatures are lacking.] 
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FIG. 4. Fit to the elastic energy loss coefficient, Q 1 (T), 
for samples of several donwgs, x ~ 0.02 (squares), 0.03 (cir- 
cles), and 0.06 (triangles) £3 All fits use the same distribution 

Of Waneiast- 

All three distributions in Fig. [l] appear to be universal, 
in that the same distribution holds over a wide doping 
range. However, different properties appear to be char- 
acterized by different distribution functions. We suggest 
that the different experiments may sample different at- 
tributes of the clusters. Thus, as discussed above, the 
susceptibility is proportional to the number of clusters 
with a particular activation energy. Since anelastic re- 
laxation responds to the change in elastic energy when 
one type of domain grows at the expense of another do- 
main with less favorable orientation, w an eiast should be 
sensitive to the domain walls, namely the perimeters of 
the domains; similarly wnqr(E) should measure the to- 
tal number of spins, or the total areas of all the domains 
with a given activation energy E. 

The nearly constant Q _1 at low T implies a nearly 
constant distribution w an eiast at low E, Fig. [j]; this can 
be understood as follows: since 
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so when Q~ x is roughly constant, so is w ane i as t(E(T)). 
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FIG. 5. Fit to NQR data Tf 1 for several dopings, x ~ 0.02 
(squares), 0.03 (circles), and 0.06 (triangles). The data are 
from Refs. |ll] and All fits use the same distribution of 
w n qr. 



III. DISCUSSION 

The distributions, Fig. [l], and the spin glass transition 
temperatures T g are the principle results of this work. As 
discussed above, w susc measures the cluster number dis- 
tribution, w ane i ast the perimeter distribution, and wnqr 
the area distribution. Figure [l] suggests that the perime- 
ters and areas follow the same distribution at high en- 
ergies, but that the area distribution starts to fall off 
sooner. 

An important finding is the universality of the dis- 
tributions: the same distributions appear to hold over 
a wide doping range. Earlier La NMR measurements 
found a similar broad distribution (fit to a Gaussian)li3, 
and found a universal weight distribution over an even 
wider doping range, including spin-frcczing^ncar-optimal 
superconducting, and Li-doped samples.ta Indeed, the 
present NQR distribution can be fit to a Gaussian, 
dashed line in Fig. [I], 
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with z = E/E m , wo — 0.367, za = 0.190, and zg — 
0.665. 

It is possible to use the distribution functions to under- 
stand how the spin slowing down evolves with tempera- 
ture, since (at any given frequency) the processes with 
the highest values of E(T) will freeze out first. Thus 
the distribution w(E) can be alternatively considered as 
a distribution w(T) of the number of processes freezing 
out at temperature T, with T ~ E/k&. From the neu- 
tron scattering, at T on (w) diagonal stripes start to order, 
with growing coherence length, but the coherence length 
cuts off at a fairly small value without diverging. At the 
same time the number of free clusters is decreasing (from 
w S usc)- There is a large change in the number of clusters 



(wsusc), but only a small decrease in the total number 
of fluctuating spins (w area = wnqr), so presumably the 
smallest, most strongly pinned clusters (largest E(T)) 
are freezing out. This interpretation that the clusters 
are small is strengthened by the fact that they make ap- 
proximately the same contribution to w per i m = w an eiast 
and Warea = wnqr (i.e., all spins are on the surface 
of the cluster). The fact that w ar ea is finite suggests 
that the spins are freezing out rather than clustering. At 
some point (perhaps when the correlation length stops 
increasing) the contribution to w area drops sharply, but 
the cluster walls can still move. Note that even though 
Warea is shrinking, there are still many freely moving 
clusters (w SU sc remains constant), suggesting that some- 
thing like a 'percolating backbone' has frozen out, and 
at lower temperatures the backbone grows by capturing 
additional clusters. 

The observed curvature in the plot of ln{w) vs. 1/T 
strongly suggests the existence of a well-defined finite 
temperature CSG transition at fixed T g . Additional mea- 
surements, especially with well-defined low frequencies, 
would help pin down the value of T g (x), but should not 
change the overall picture. The interpretation is com- 
plicated by the presence of a broad distribution of acti- 
vation energies, which lead to slowing down over a wide 
frequency range. In particular, we have shown that an 
apparent scaling of the order parameter q can arise from 
a particular form of w{E), with sharp upper cutoff E m . 
However, scalmg behavior is also found in an applied 
magnetic field,l3 and this is harder to explain. The dop- 
ing dependence of E rn (insert, Fig. |^) is very suggestive 
of pinning effects. Note that E m appears to diverge as 
x — > xo ~ 0.01, but this is cut off by the Neel transition 
near x — 0.02. We suggest that the observed behavior 
is characteristic of stripe freezing in the presence of de- 
fect pinning. In the absence of pinning, there could be 
a long range stripe ordering phase transition, but only 
at T — > due to fluctuations of the walls in two dimen- 
sions. In this case, the correlation length would diverge 
as T — > 0. Impurities will tend to pin the charge stripes, 
but the pinning will vary in space, depending on the dis- 
tribution of the impurities, and hence leading to a wide 
distribution of activation energies. Certainly, the large 
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(Fig. [ 

- J/4, 



) is sugges- 
with J the 



tive of an electronic process 
exchange constant. 

In attempting to refine the experimental determina- 
tion of T g (x), it must be noted that T on may be a strong 
function of the sensitivity of the probe to fluctuations. 
For instance, whereas the La NQR shows a T on in the 
range of ~ 10-20 K, Cu NMR is wiped out at a much, 
higher temperature, ~ 50 K. It has been suggesteal3'E£l 
that this wipeout is caused by the same CSG fluctua- 
tions, which are in the Cu02 plane, and hence have a 
much greater effect on Cu than on La (for an alternative 
interpretation, see Ref. |l9| ). A remaining puzzle is that 
the wipeout should start to recover below T g jEB whereas 



no recovery is found for the Cu signal down to 350 mK.I 
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IV. CONCLUSION 

The signatures of spin freezing in La2_ x Sr a; Cu04 from 
neutron spectroscopy, NQR relaxation, anelastic relax- 
ation and magnetic susceptibility (spanning ~ 13 orders 
of magnitude of frequency) have been reconsidered in 
terms of a distribution w (E) of activation energies for 
the spin dynamics. Each experiment probes a different 
aspect of the spin clusters, like the total number of spins 
in a cluster, the spins along its perimeter or the number 
of clusters with an unpaired spin, and therefore requires a 
different weight function for the distribution of activation 
energies. It is found, however, that for each process the 
same w (E / ' E m (x)) can explain all the data, independent 
of doping, with the same E m (x) for all the experiments. 
In addition, E m (x) scales linearly with the temperature 
T g (x) at which the effective spin fluctuation rate deduced 
from all the experiments freezes toward the cluster-spin 
glass state. 
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